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Abstract 



A Monte Carlo simulation study of the vacancy-assisted domain-growth in 
asymmetric binary alloys is presented. The system is modeled using a three- 
state ABV Hamiltonian which includes an asymmetry term, not considered 
in previous works. Our simulated system is a stoichiometric two-dimensional 
binary alloy with a single vacancy which evolves according to the vacancy- 
atom exchange mechanism. We obtain that, compared to the symmetric case, 
the ordering process slows down dramatically. Concerning the asymptotic 
behavior it is algebraic and characterized by the Allen-Calm growth exponent 
x = 1/2. The late stages of the evolution are preceded by a transient regime 
strongly affected by both the temperature and the degree of asymmetry of 
the alloy. The results are discussed and compared to those obtained for the 
symmetric case. 
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I. INTRODUCTION 



The study of domain-growth is a prototype problem in out-of-equilibrium statistical 
mechanics. Besides its fundamental interest, it has many technological implications in areas 
like metallurgy, semiconductors, surface physics, ...etc. Despite during the last 15-20 years 
domain-growth has been the subject of a deep and continuous investigation |1].@,[|, the 
models commonly used are still far from real materials. A typical experimental situation 
occurs when a binary alloy undergoing an order-disorder transition is rapidly cooled from 
the high-temperature disordered-phase to the low-temperature ordered-phase. The growth of 
the ordered domains, subsequent to the quench, is a complex phenomenon in which different 
factors may play an important role: details of the atomic diffusion, topological defects, 
temperature, degeneration of the ordered phase, interface structure, impurities, etc... It is 
important to realize that the information relative to the influence of each single factor is 
not always available from the experiments and therefore alternative reliable methods such 
as numerical simulations may be very useful. 

Many Monte Carlo simulation studies based on the Ising-like lattice models assume a 
very simplified picture of the alloy: i.e., perfect, without defects and with an ordering 
mechanism based on the neighboring atom-atom exchange. The exchanges proposed are 
accepted or rejected according to the Metropolis algorithm. The results obtained are 

in agreement with the phenomenological Allen-Cahn curvature-driven theory [|7| for the time 
evolution of the mean-size of the ordered domains, i.e. R(t) ~ t x , with x = 1/2. Moreover, 
evidence has accumulated that the ordering process exhibits dynamical scaling and that the 
late-time growth behavior is subjected to a high degree of universality |§. Concerning the 
experiments, only a small number have been designed with the aim to check the validity 
of the growth-law and one may say that, in general, they are compatible with the growth 
exponent x = 1/2, at least, for quenching temperatures below but close to the equilibrium 
ordering temperature. In reference || the authors report a very complete summary of the 
experimental data available in the literature. 

Concerning binary alloys, atomic diffusion is the underlying microscopic mechanism in- 
herent to any ordering process. Moreover, it is accepted that the atomic diffusion occurs via 
vacancies. In this sense, quite recently, the possibility for the more realistic vacancy-assisted 
domain-growth has been investigated [|T^,|n|,|T^] . The model description is based on the 
three-state ABV lattice model |J and the microscopic dynamical mechanism accounts for 
the atom-vacancy exchanges only. The most important result is the prediction for growth 
exponents greater than the standard Allen-Cahn value. 

In this paper we make another step towards the study of real alloys. In most theoretical 
studies it has been assumed that the alloys are symmetric, i.e. the Hamiltonians used 
are invariant under the exchange of the two species, A and B, forming the alloy. In the 
framework of pure Ising models with pair interactions the asymmetry term appear as a field 
term (chemical potential) which is usually neglected. This turns out to be irrelevant for the 
models without vacancies since the conservation law for both A and B particles makes the 
field term to be constant. Nevertheless, for models including vacancies, such asymmetry 
terms should be considered. 

The motivation for this study comes from the fact that in many binary alloys the va- 
cancies exhibit a tendency to locate preferentially in one of the possible different equivalent 
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sublattices (into which the ordered structure can be decomposed) [[L3|]. Our main goal is to 
modelize such behavior by taking into account the corresponding asymmetry term in the 
ABV model, and perform Monte Carlo simulations of the vacancy assisted domain growth. 
Although the present simulations are restricted to a two-dimensional lattice, we expect that 
the general conclusions will apply to real alloys. 

The results obtained may be summarized as follows: The main effect of the asymmetry 
parameter is to decrease dramatically the speed of the ordering process. The specific inter- 
action between vacancies and interphases is repulsive and, therefore, the vacancy migration 
from the ordered bulk to the interphases, which is required for the domains to grow, needs to 
be thermally activated. The associated energy barrier gives rise to a temperature-dependent 
transient previous to the true asymptotic behavior. Moreover, the process of ordering inside 
the domains lasts until long times. Nevertheless, we obtain that the long-time domain- 
growth behavior is algebraic, and compatible with the Allen-Calm theory with an exponent 
x = 1/2 at all the temperatures and degrees of asymmetry studied. This result is different 
from previous studies of vacancy-assisted dynamics of domain-growth in symmetric binary 
alloys where the standard x — 1/2 value is found only at temperatures T —>■ T c while, at 
moderate and low temperatures, the exponent is definitively much larger than 1/2 and tends 
to 1 for T -> 0. 

The remaining of the paper is organized as follows. Next, in section II, we briefly sum- 
marize the ABV model Hamiltonian and discuss the behavior of the vacancies in terms of 
the model parameters. In the same section we also define the model dynamics. In section 
III we present the details of the Monte Carlo simulations and provide the definition of the 
relevant quantities computed. Results are presented and discussed in section IV. Finally, in 
Section V we conclude. 



II. THE MODEL 

A. The ABV Hamiltonian 

The ordering configurations of an AB alloy with a constant concentration of vacancies, 
are described in terms of three state spin-1 variables Si = 1, —1,0, defined on each lattice 
site i = l, 2....N. Si = 1 means that site % is occupied by an A atom, Si — — 1 by a B atom 
and Si = by a vacancy. We shall restrict to 2d-square lattices with a constant number 
of particles (Na,Nb) and vacancies {Ny). Considering only nearest neighbor (n.n.) pair 
interactions, the Hamiltonian for the ABV model can be written as PJHf: 

n = JY,s*s 3 + ky,s?s] + + s^) (1) 

(*?') (ij) (ij) 

where the sums extend over all n.n. pairs. It contains three independent parameters: J, K, 
and L. The parameter J determines the ordering of the system. We take J > in order 
to ensure that the ground state will be antiferromagnetically ordered and formed by two 
alternating sublattices (A-sublattice and B-sublattice). In what follows we shall take J as 
the unit of energy. The parameter L* = L/J accounts for the energy difference between 
AA and BB bonds, i.e.. it controls the asymmetry of the alloy. Finally, the parameter 
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K* = Kj J controls the energy of the bonds involving vacancies, and therefore the vacancy- 
vacancy specific interaction. In Table | we have summarized the different bond energies and 
its corresponding excess-energy with respect to the ordered AB bond. 

The behavior of the vacancies during the ordering process depends on both parameters 
K* and L*, and may be understood from the following considerations: 

1. For given values of Na, Nb{= Na) and Ny, the asymmetry parameter L* controls the 
tendency for the vacancies to locate preferentially in one of the two sublattices. If, for 
instance, L* > 1 (L > J) the energy of the BB bond is lower than the energy of the 
AB bond. The vacancies, in particular those initially located in the A-sublattice, will 
show a tendency to migrate to the B-sublattice. The same argument holds if L* < — 1 
for the A-sublattice. 

2. Furthermore for \L*\ > 1, the tendency for the vacancies to concentrate in a given sub- 
lattice, determines the specific interaction between vacancies and antiphase boundaries 
(APB). The APB's are interphases separating thermodynamically equivalent ordered 
domains with the same absolute value of the corresponding order parameter. These 
APB's are always present during the domain-growth process and are formed by a se- 
quence of AA and BB bonds. For values of \L*\ < 1, the AA and BB bonds are 
unfavorable with respect to the AB bond, and therefore the vacancies tend to concen- 
trate at the APB's. But if, for instance, L* > 1 the BB bond is energetically the most 
favorable and the vacancies, located preferentially in the B-sublattice, will not show 
any natural tendency to go to the APB's (notice that the AA unfavorable bonds in 
the APB's can only be broken by vacancies located in the A sublattice). In short, for 
the asymmetric case (\L*\ > 1), the vacancies will prefer to stay in the ordered bulk 
rather than to concentrate at the APB's, and the corresponding energy difference is 
proportional to (\L*\ — 1). 

3. Finally, we notice that the behavior described above may be modified by the specific 
interaction between vacancies, controlled by the parameter K*. A simple analysis 
shows that this interaction is attractive for K* < 1 and repulsive for K* > 1. Recently, 
the effect of K*, in the simplest case L* = 0, has been extensively studied |15| , [l6| , [r7| , [18| . 
Here we shall concentrate on the effect of L* only and leave for a further investigation 
the interplay between K* and L* that may be very subtle. 

From the above considerations one finds six different regions in the space of model pa- 
rameters, as it is illustrated in Fig.|I|. We notice that the different behaviors of the vacancies 
are separated by the lines \L*\ = 1 and K* = 1. 



B. Model dynamics 

The non-equilibrium properties of model (|I|) may be studied using Monte Carlo com- 



puter simulation techniques [19]. This requires the implementation of a given microscopic 



ordering mechanism, compatible with the conservation laws in effect. In our case this means 
to preserve the number of particles of each kind as well as the number of vacancies. We 
have chosen the vacancy mechanism: a single vacancy is considered and only atom-vacancy 
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exchanges are allowed. We call this a vacancy jump. Proceeding as it is usual, the vacancy 
jumps are proposed to n.n. and next nearest neighbors (n.n.n) sites, with equal probabil- 
ity. This is very convenient when performing Monte Carlo simulations in order to prevent 
trapping phenomena P,|T0|1. The attempts are accepted or rejected following the Metropolis 
algorithm which only takes into account the energy difference between the initial and final 
states in a given vacancy jump, that is: P(A7i) = min{l, exp( ~n^ )}. The unit of time, the 
Monte Carlo step (MCS), is defined as N attempts of single vacancy jumps. 

In the present study the word "barrier" is used to denote the extra energy needed for 
the vacancy to migrate from the bulk to the antiphase boundaries. This barrier is intrinsic 
to model (|I]) and should not be confused with eventual energy barriers associated to the 
details of the vacancy path in a given jump between two neighboring sites. There exist in 
the literature Monte Carlo algorithms which consider the existence of these last barriers 
f20|| . Nevertheless, for the sake of clarity, in this first study of domain growth in asymmetric 
binary alloys, we will not consider them. 



III. SIMULATION DETAILS 

Our simulated system is a (nearly) stoichiometric AB alloy containing a single vacancy. 
Notice that, in this case (Ny = 1, Nb = Na — 1), the parameter K* is irrelevant and 
therefore may be neglected in model (H). The particles (N — 1) are sitting on a square lattice 
of size i x i = N, subjected to periodic boundary conditions. The main results presented 
correspond to L* = 2 and I = 600 (N = 3.6 • 10 5 ). Some additional results for other values 
of the asymmetry parameter (L* = 0, L* = 1 and L* = 3) are also shown. Starting from a 
completely disordered configuration, the ordering process is studied at different quenching 
temperatures (T* = ksT/J = 1.0, 1.25 and 1.5) below the critical ordering temperature 
T* ~ 2.27. The simulations are extended up to ~ 4 • 10 4 MCS and the results are averaged 
over ~ 250 independent runs. 

From the information given above it follows that the concentration of vacancies in our 
simulated system is ~ 10~ 6 . This is reasonable since typical concentrations of vacancies may 
range from ~ 10~ 3 to ~ 10~ 10 , strongly depending on temperature | 2~lf . For the case with 
L* = 0, it has been shown that the results obtained in vacancy-assisted dynamics do 
not change if one increases the number of vacancies (and the system size) maintaining its 
concentration constant. This is because of the small vacancy concentration we are dealing 
with which makes the vacancy- vacancy interaction term in eq. ([TJ) to be negligible in front of 
the other contributions. We expect that the same will apply in the present non-symmetric 
study and, therefore, focus on the simplest case Ny = 1. 

From the simulations we have measured the following quantities: 

1. Long range order parameter m. It is defined as the absolute value of the antifer- 
romagnetic order parameter on the system. 

1 N 

m = — SiSign(i) (2) 

iV i=l 

where the function sign(i) takes alternating ±1 values on the lattice in a chessboard 
way. 
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In order to obtain information concerning the path of the vacancy during the ordering 
process it is convenient to define the two following local order parameters: 



2. Local order parameter around the vacancy m v . It is defined as the absolute 
value of the antiferromagnetic order parameter in a travelling (5 x 5) cell centered at 
the vacancy position. In order to reduce the numerical uncertainties partial averages 
over consecutive MCS have been performed. In the case of L* = 2 each value is the 
average over 10 (consecutive) MCS while for L* = 3 the same has been done over 
100 MCS. Furthermore, averages over many (30-1000) independent runs have been 
performed in order to minimize the dependence with the initial conditions. 

3. Local order parameter m( 5x5 ). It is defined as the absolute value of the antiferro- 
magnetic order parameter in a (5 x 5) cell. The computation is done by dividing the 
original lattice into cells of size (5 x 5) and then averaging over all of them. 

Domain-growth can be monitored by measuring other quantities. We shall study the 
excess-energy and the width of the structure factor. 

4. Excess internal energy AE(t). It is defined as the excess internal energy per 
particle: 

AE(t) = E(t)- <H> T (3) 



where E(t) is the energy of the system (given by eq. (jl])) at time t after the quench 
and < Ti. >t is the equilibrium energy of the system at the temperature T. At 
each temperature, this equilibrium energy has been obtained on a system of i = 200, 
starting from a completely ordered configuration and then using the standard atom- 
atom exchange mechanism. The final value of < Ji >t is the average over the last 
1.8 • 10 4 MCS of a single evolution of 2 ■ 10 4 MCS long. 

5. Structure factor. The structure factor is defined as the Fourier transform of the 
correlation function. It is written as: 

2 

(4) 



S(k) 



1 

N 



E^exp 



where k is a reciprocal space vector, a is the lattice spacing, and fl is the position 
vector of site %. We have computed the profiles along the (10) and (11) (and equivalent) 
directions, around the superstructure peak at k = (\\)- 

6. Structure factor width u{t). It is defined as the square root of the second moment 
of the structure factor: 

Q m ax / Qmax \ 

E <? 2 %) / E %) (5) 

q=0 I q=0 J 
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where q = \k — (||)| is the distance to the superstructure peak. The sum is performed 
over the q values in the corresponding direction and extends up to q max = 3a(t). Note 
that a {t) and q max have been obtained in a self-consistent way. By this method we avoid 
the problems associated to the fitting of a profile function and to the determination 
of the background which, in this problem, cannot be neglected and shows a time 
evolution. 



Then, provided scaling holds, 

AE(t) ~ m(t) _1 ~ a (t) ~ FC l {t) (6) 

where R(t) is the mean size of the ordered domains. Finally, in order to reveal the existence 
of transient regimes during the growth, it is convenient to introduce a definition for the 
growth exponent, not affected by a priori assumptions: 

7. Effective growth exponent x e (t). It is defined as the logarithmic derivative of the 
excess-energy according to: 

= -rfbgtAEW) = -bgtAfl^l/AEW) 
v 1 dlog(t) v-*o logz/ w 



where, in the computation, instead of the limit, we have taken a small value of v 
[y = 2) and fitted a power-law to all intermediate data contained in the range (t, ut). 



IV. MONTE CARLO RESULTS AND DISCUSSION 



In this section we present the numerical results obtained from the Monte Carlo simu- 
lations. Figure ^| shows snapshots of the configurations obtained at T* = 1, for L* = 
and L* = 2 at four different selected times. Disordered sites are shown in black. The 
snapshots for the two values of L* correspond to simulation times for which the system has 
the same excess internal energy. The overall picture of the domains is qualitatively similar. 
Nevertheless, for L* = 2 one observes that the amount of disorder in the bulk is larger. 

Prior to the study of the growth-law itself, it is interesting to analize the behavior of 
the vacancy during the ordering process. More precisely, it is important to understand the 
characteristics of the path followed by the vacancy during the vacancy-assisted ordering 
process which takes place as a response to the quench, performed at T* . This is done by 
monitoring the behavior of m v vs. m( 5x5 ) as it is illustrated in Fig. |3| for a system of linear 
size i = 200, temperature T* = 1.0 and four different values of the asymmetry parameter 
L*. The straight line denotes a random walk eventually followed by the vacancy in case 
the path is uncorrelated with the state of order in the system. As it can be seen in Fig. 
|!| for L* = 2 and L* = 3, the order around the vacancy is clearly higher than the average 
local order in the system (at all times m v > 777.(5x5))- This is indicative of the tendency of 
the vacancy to stay preferentially in the ordered regions (inside the bulk). This behavior is 
opposite to the one found for L* = |l0|Jlll| , where the vacancy clearly exhibits a natural 
tendency to concentrate at the interphases (APB's). Notice that, after an initial transient, 
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the intermediate case L* = 1 corresponds to the uncorrelated random walk. Moreover, it is 
clear that the process rapidly slows down as L* increases (especially for \L*\ > 1). This can 
be realized from the number of MCS (indicated by an arrow) needed to achieve the same 
amount of order in the system (for instance a value of m( 5x5 ) = 0.50). Notice also, that 
in all cases, at short times, m v > m(5 X 5)- This is because the ordering mechanism involves 
the vacancy jumps exclusively and therefore, at short times, the system can only be locally 
ordered around the vacancy. The above scenario is in full agreement with the discussion 
given in section II i.e., for \L*\ > 1 the specific interaction between the vacancy and the 
APB is repulsive. In what follows we shall mainly concentrate in the asymmetric case with 
L* = 2. 

We first study the time evolution of the growing domains by following the behavior of 
AE(t). This is shown, in a log-log plot, in Fig. ^for three different quenching temperatures 
T* = 1.0, 1.25, 1.50 and £ = 600. In order to verify that finite size effects do not mask the 
long time behavior, we have also measured the corresponding time evolution of the order 
parameter m(t), which is shown in Fig. |5|, for the same three temperatures. Note that 
simulations have been run up to a time (~ 4 • 10 4 MCS) for which the order parameter is 
m ~ 0.1. This is enough to ensure that, in the long time stages studied, the system is still 
in the regime of competing domains. 

Results in Figures |] and [5] indicate that the asymptotic regime of AE(t) and m(t) 



22] is, for the three temperatures, compatible with the algebraic Allen-Cahn growth-law, 



denoted in the figures by a dotted straight line. We have estimated the growth-exponents 
of the different curves by least-square fitting a power-law to the last ~ 25000 MCS. They 
are indicated in the insets. All the values are compatible with the Allen-Cahn exponent 
(x = 0.50 ± 0.06). Nevertheless, notice that the values of the exponent obtained from 
AE(t) show a slight tendency to increase as one decreases the quenching temperature. This 
has to do with the temperature-dependent transient regime, not revealed by the behavior 
of m(t), that may affect the asymptotic behavior up to very long times. Moreover, from 
our simulations we have observed that, inside the domains, the equilibrium ordering is 
established asymptotically and is not complete until long times (see Fig. ^|). From all these 
considerations it follows that, in our case, AE(t) has to be taken as an approximate measure 
of the domain size evolution and that the corresponding exponents are only orientative of 
the tendency of the system. We now turn back to the transient regime. It is related to the 
existence of thermally activated processes which hinder the growth of the ordered domains. 
Indeed, once the domain structure is formed, the domains grow only when the vacancy is 
located at the APB's. Since for \L*\ > 1 the vacancy prefers to stay in the ordered bulk (as 
it was discussed in section II), the migration to the APB's needs to be thermally activated. 
The corresponding energy barrier that the vacancy has to climb can be easily evaluated 
and one obtains e* = 4(|L*| — 1) = 4. This barrier verifies e* > T*, even for the highest 
temperature studied. 

A more quantitative analysis of the whole evolution can be extracted from the compu- 
tation of the effective exponent obtained from the logarithmic derivative of AE(t), as it 
was defined in the previous section. Results are shown in Fig. |[ The general trends are, 
for the three temperatures studied, the same. The effective exponent first increases from 
very small values (at short times) up to a maximum (over x = 1/2) at a time to and next 
decreases to reach, asymptotically the Allen-Cahn value. On the other hand, the position of 
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the maximum (to) and the value of the maximum depend on temperature: both are larger 
as lower the temperature. We have evaluated the value of to for the different temperatures 
and get t = 1 • 10 3 , 0.59 • 10 3 and 0.64 • 10 3 MCS for T* = 1.0, 1.25 and 1.50 respectively. 
Moreover, the value of t depends on L*. In fact, we have performed simulations for L* = 3 
(not shown in the paper) at the same temperatures, and get to — 4.1 • 10 4 , 1 ■ 10 4 and 0.58 ■ 10 4 
MCS respectively. The value of to is, for L* = 3 much larger than for L* = 2, providing an 
idea about the increasing times needed to reach the asymptotic regime when increasing L*. 
We expect that to is proportional to the passing time of the barrier, 

to ~ e e */ T * (8) 

We can check that e* oc L* — 1 by representing to vs. the factor (L* — 1)/T* in a linear- 
log plot. This is shown in Fig. |7| for different values of L*. The square corresponds to 
L* — 1, the diamonds to L* = 2 and the dots to L* = 3. The dashed line shows an 
exponential fit, which renders e* = 3.7(L* — 1), in agreement with the theoretical estimation 
given previously, e* = 4(|L*| — 1). The deviation of the points corresponding to the highest 
temperature (T* = 1.5) may be attributed to the existence of higher order corrections in 
Eqn. ([5p due to entropic effects. 

We now present the data obtained from the study of the structure factor. In Fig. || we 
show the time evolution of the width of the superstructure peak, a(t), in the two relevant 
directions and for the same three temperatures and model parameters as in Figures |] and |5|. 
The long-time behavior is, for the three temperatures, algebraic (a ~ t~ x ) and characterized 
by an exponent x = 0.50 ± 0.05. The inset contains the values of the growth exponent 
obtained by fitting the algebraic law a(t) ~ t~ x to the data at each temperature. 

Moreover we have studied the existence of dynamical scaling. As an example, in Figure 
|9]we show the structure factor profiles at T* = 1.0 scaled with the corresponding a(t) along 
the two directions. The overlap of the data is very satisfactory. Moreover the shape of 
the scaling function can be compared with the universal scaling function corresponding to a 
non-conserved order parameter system p3[, which is also plotted with a continuous line. The 



agreement is excellent for large and intermediate values of q/a. For small values q/a < 1 
the behavior is slightly different in both directions. This may be associated to the existence 
of anisotropic correlations in the domain shape. This will be analyzed in detail in a future 
work. 



V. SUMMARY AND CONCLUSIONS 

The results presented in the previous section show that during the asymptotic regime 
of vacancy-assisted domain-growth A£'(t) _1 ~ m(t) ~ alt)^ 1 . In spite of some details 
observed for small values of q, one may conclude that dynamical scaling holds and that 
the long-time behavior can be characterized by an unique length which evolves with time 
according to R(t) ~ t 1 ^ 2 . This asymptotic regime is preceded by a long transient that 
depends on both the temperature and the asymmetry character of the alloy. In the ABV 
model the asymmetry is controlled by the parameter L* which in turn determines the energy 
barrier e* = 4(|L*| — 1) that the vacancy has to climb in its migration from the bulk to the 
interphases. The characteristics of the transient regime are determined by the ratio e*/T*. 
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This behavior is definitively different from the one found in symmetric binary alloys with 
L* = 0. The origin of this difference lies in the fact that for symmetric binary alloys the 
interaction between the vacancies and the interphases is attractive. The first consequence 
of this is to speed up the global process. The growth exponent was found to be temperature 
dependent || and clearly larger that 1/2 at low and moderate temperatures. Only at 
T* — > T* the Allen-Cahn exponent is recovered. 

The fast growth predicted in the symmetric case has never been found experimentally. 
The authors ||10| , |rTf argued that the temperatures at which experiments are usually per- 
formed are too high. In the light of the present results it seems more reliable that it has to 
do with the intrinsic asymmetric character of real alloys. Actually the asymmetry parameter 



L* can be estimated from ab-initio calculations. We have used the WIEN97 code 24 to 



evaluate L* for the /3-CuZn alloy and have obtained L* = 3. p5 
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TABLES 

TABLE I. Absolute bond energies and excess energies with respect to the A — B bond for the 
ABV model as a function of the parameters J, K, and L. 

Bond Energy Excess energy 

A-B -J + K 

A- A J + K + 2L 2J + 2L 

B-B J + K-2L 2 J - 2L 

A-V J -K 

B-V J -K 

V-V J-K 
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FIG. 1. Space of the parameters K* and L* divided into the regions where different dynam- 
ics are expected. The solid line K* = 1 separates the region of vacancy attraction (K* < 1) 
from the region of vacancy repulsion (K* > 1). The solid lines L* = ±1 separates the region 
where the vacancies tend to precipitate at the antiphase boundaries (\L*\ < 1) from the region of 
vacancy-antiphase boundary repulsion (\L*\ > 1). 
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FIG. 2. Sequence of snapshots of the evolving domain structure for a system with £ 
= 1.0 and two values of L*. Black regions represent disordered sites 
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FIG. 3. Evolution of the local order parameter defined around the vacancy m v vs. the local 
order parameter ms X 5 for a system with £ = 200, T* = 1.0 and different values of L*, as indicated. 
The dotted straight line shows the behavior corresponding to a random walk of the vacancy. The 
arrows indicate, for each case, the number of MCS necessary to reach a value m5 X 5 = 0.5. 
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FIG. 4. Log-log plot of the time evolution of the excess-energy AE(t) for a system with L* = 2, 
I = 600 and three different temperatures. The dotted straight line indicates the Allen-Cahn 
behavior AE ~ t~ l l 2 . The inset shows the fitted power-law exponents to the last 25000 MCS. 
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FIG. 5. Log-log plot of the time evolution of the long-range order parameter m, for the same 
temperatures than in Fig. 4. The dotted straight line indicates the Allen-Calm behavior m ~ 
The inset shows the fitted power-law exponents to the last 25000 MCS. 
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FIG. 6. Effective growth exponent as a function of time obtained from the logarithmic derivative 
of the behavior of AE(t), corresponding to the same cases as in Fig. 4. Estimations of to an d 
typical error bars are shown. 
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FIG. 7. to as a function of (L* — 1)/T*. The dashed line shows the fitted exponential behavior 
with a slope e*/(L* - 1) = 3.7. 
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FIG. 8. Evolution of the structure factor width a(t), for a system with L* = 2 and I = 600, 
along the two directions (10) (continuous lines) and (11) (discontinuous lines). Data corresponds 
to three different temperatures, as indicated. The dotted straight line indicates the Allen-Cahn 
behavior a(t) ~ t" 1 / 2 . The inset shows the fitted power-law exponents to the last 25000 MCS. 
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FIG. 9. Scaling of the structure factor profiles for a system with L* = 2, I = 600 and T* = 1.0, 
along the two directions (10) and (11). The continuous line corresponds to the universal scaling 
function proposed by Ohta et al. (Ref. 23). 
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